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We investigate the polarization modulation properties of a variable-delay 
polarization modulator (VPM). The VPM modulates polarization via a 
variable separation between a polarizing grid and a parallel mirror. We find 
that in the limit where the wavelength is much larger than the diameter 
of the metal wires that comprise the grid, the phase delay derived from 
the geometric separation between the mirror and the grid is sufficient to 
characterize the device. However, outside of this range, additional parameters 
describing the polarizing grid geometry must be included to fully characterize 
the modulator response. In this paper, we report test results of a VPM at 
wavelengths of 350 //in and 3 mm. Electromagnetic simulations of wire grid 
polarizers were performed and are summarized using a simple circuit model 
that incorporates the loss and polarization properties of the device. © 2011 
Optical Society of America 

OCIS codes: 350.1270. 120.5410, 230.4110, 240.5440, 050.6624. 
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1. Introduction 


Astronomical polarimetry in the far-infrared through microwave portion of the electromag- 
netic spectrum is a useful tool for probing the physics of interstellar dust, investigating the 
role of magnetic fields in star formation, and characterizing the radiation from the early uni- 
verse. In each of these cases, polarization modulation is an important element of instrument 
design because it enables a precision measurement by encoding the polarization information 
and thus separating it from the typically larger unpolarized background signal. 

Variable-delay polarization modulators (VPMs) operate through the introduction of a 
controlled, variable phase delay between two orthogonal linear polarizations. This is accom- 
plished in recent work by placing a polarizing grid in front of a mirror and varying the 
separation [1,2], although other related architectures have been used for polarization mod- 
ulation [3,4]. For a single frequency, the polarization transfer function of the VPM is a 
sinusoidal function of the phase delay between the two orthogonal polarizations. The key 
to understanding the polarization transfer function is to determine how the phase delay is 
related to the grid-mirror separation, since the latter is the quantity that is directly measur- 
able in an instrumental setup. In practice, the polarization response of a specific modulator 
can be measured to produce a template for use in reducing astrophysical data. As such, the 
work presented here is not required for demodulating a time-ordered data set. The utility 
of the model described here is that it allows calculation of the achievable optical response 
and provides both guidance regarding the modulators sensitivity to design parameters and 
validation metrics for the observed instrumental performance. In addition to VPMs, grid- 
mirror systems have been implemented in a variety of other applications [5-11]. The analysis 
presented here is applicable to these systems as well. 

In this paper, the polarization transfer function for a single modulator from the 
Hertz/VPM instrument is measured at 350 /an and 3 mm. The observed response is well- 
characterized by a transmission line model. This work represents a refinement of previous 
results [2] where it is qualitatively shown that VPMs modulate polarization with reason- 
able efficiency and without introduction of major artifacts. Here, we present new lab tests 
of Hertz/VPM collected after improving the optical alignment. We demonstrate quantita- 
tive agreement between the measured polarization transfer function and a transmission line 
model. Also, we show that the overall modulation efficiency of Hertz- VPM is within ~2% 
of the expected value. The organization of the paper is as follows: In section 2, we briefly 
review the principle of operation of the VPM. In section 3 we describe the circuit model for 
the VPMs. In section 4, we discuss the current state-of-the-art for analytical models used 
to calculate the circuit parameters and include numerical simulations to supplement then- 
limitations. In sections 5 and 6, this model is applied to measurements of the VPM. We 
discuss the results in section 7. 
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2. The Variable-delay Polarization Modulator (VPM) 


Input Port Output Port 



Fig. 1. The VPM consists of a polarizing grid placed in front of and parallel 
to a mirror. The polarization parallel to the grid wires is reflected by the grid. 
The orthogonal linear polarization passes through the grid and is reflected off 
the mirror. The two components are recombined at the output port with an 
introduced phase that is dependent upon the grid-mirror separation, d. 


The VPM consists of a polarizing grid positioned parallel to and in front of a mirror. A 
change in grid-mirror separation corresponds to a change in introduced phase between two 
linear orthogonal polarizations. The device is show in Figure 1. Choosing coordinates such 
that Stokes Q gives the difference between the polarization states parallel and perpendicular 
to the VPM wires, the polarization transfer function can be expressed as 

U' — U cos 5 + V sin 5. (1) 


Here, U and U' are the input and output Stokes U parameter, and V is the input circular 
polarization. Here, 8 is the electrical phase delay between the polarized components trans- 
mitted by and reflected by the grid. This is the phase of interest when using the VPM as a 
modulator. 

In the limit for which the wavelength is much larger than the length scales that characterize 
the local grid geometry, the VPM phase is proportional to the path difference between the 
two polarizations, 


^oo 


And 

X 


cos 9. 


( 2 ) 


where 0 is the incident angle, as shown in Figure 1. Here, the oc subscript indicates that this 
relationship holds in the long wavelength limit. In many polarizing grid applications the size 
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of the wires relative to wavelength is large enough to produce measurable deviations from 
the phase relation given in Equation 2. 

We assume throughout this paper that the grid wires are uniform, infinite in extent, and lie 
in a plane. Diffraction grating lobes are negligible when the wire radius, a, and grid constant 
or center-to-center wire pitch, g, are much smaller than the wavelength of the incident 
radiation. In this limit, the structure has a homogeneous response for each polarization 
mode assuming plane-wave illumination. 


3. A VPM Model 


An “ideal” polarizing grid perfectly separates orthogonal linear polarization components 
with a response that is independent of frequency. That is, the grid perfectly transmits the 
polarization component having the electric field perpendicular to the grid wires and perfectly 
reflects the component having the electric field parallel to the grid wires. The former corre- 
sponds to zero impedance contrast between free space and the grid; the latter corresponds 
to an infinite impedance contrast. For physically realizable grids, the achievable impedance 
contrast for each of the two polarizations depends on the electromagnetic properties and 
geometry of the structure. In the limit the ware radius and separation are finite compared 
to the wavelength of the incident radiation, the response becomes dependent on frequency. 
Motivated by the observation that the wire grid is a polarization- and frequency-selective 
surface, we model the structure using a transmission line approach. 

We model the polarizing grid using a circuit representation [12] of the VPM as shown in 
Figure 2. We normalize all impedances to that of free space, Zq = 377 D per square. For 
the polarizing grid, the transmission and reflection of the electric field component parallel 
to the grid wires can be represented by an inductive circuit. The transmission or A BCD 
matrix [13, 14] that represents the grid for this polarization is 
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( 3 ) 


Likewise, the grid’s effect on the polarization component having its electric field perpen- 
dicular to the grid wires can be represented by a capacitive circuit and is given by 


1 / iX c + iXcc + (2 + Re (2 Rc + iX cc )iX c 

>Xcc + (2 + R c \ 2 + Xv ( 2 + jfe) Rc + i 1 + Tcc) lXcc 

(4) 

To calculate the VPM transfer function, each of these matrices can be cascaded with the 
transmission matrices for a free space delay and a short. The resulting matrix is employed 
to calculate the reflection coefficients for each polarization, F|| and Tj_, as a function of 
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Parallel Polarization - Inductive 



Wire Grid Free-Space Delay Mirror 


Fig. 2. In the limit the wavelength is large compared to the wire pitch, the VPM 
can be modeled by two independent circuits. An inductive circuit is used for 
the polarization component having the electric field parallel to the grid wires 
(top). A capacitive circuit is used to model the polarzation component having 
the electric field perpendicular to the grid wires (bottom). 



Fig. 3. The geometry for the general polarizing grid lying in the x-y plane is 
shown. 
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frequency. The Jones matrix [15] for the VPM can then be expressed as 


J = 




( 5 ) 


The Stokes transfer function can be calculated via a density matrix. The density matrix 
may be written as a linear combination of Stokes parameters (/, Q, U, V) and can be used 
to completely characterize the polarization state of an arbitrary partially-colierent radiation 


field 



+ Q 


1 0 
0 -1 




( 6 ) 


The density matrix of the incoming light can be mapped to that of the outgoing light by 

t y - 

This model is applicable for normal incidence. In the more general case of non-normal 
incidence, it is convenient to parameterize the model in terms of the incidence angle, 9, and 
the angle between the grid wires and the plane of incidence, cp. See Figure 3. This general 
case can be calculated using the above formalism along with a procedure similar to that 
described by Goldsmith [14] where the free-space TE arid TM polarizations are projected 
onto the grid wires. The polarization separation then occurs in the basis of the wires, as 
done above. 

If the incidence angle is nonzero, there are two effects. First, the length of the grid-mirror 
separation in the transmission line is multiplied by cos#. Second, it is necessary to take into 
account the differences at the grid interface for the TM and TE polarizations. For the case 
of the transverse magnetic (TM) polarization mode, the effective impedance is altered from 
the free-space value by a factor of cos 6. This is due to the fact that the electric field that 
is projected onto the plane is smaller than the freely-propagating value. For the transverse 
electric (TE) polarization, it is the magnetic field vector that is projected into the plane and 
so the impedance differs from that of free space by a factor of 1/ cos#. 

The simplest case that one can consider with a non-zero incidence angle is that where the 
wires are either parallel or perpendicular to the plane of incidence. In this case, each of the 
TE and TM modes correspond to only one of the circuits described above. For example, if 
the wires are parallel to the plane of incidence, the TE and TM modes correspond to the 
capacitive and inductive circuit elements, respectively. For a general grid rotation angle (p, 
four circuit models must be considered utilizing all of the combinations of (TM,TE) and 
(inductive, capacitive). After the transmission matrices are found for the two transmission 
lines in Figure 2, scattering parameters can be determined. Once the scattering matrices are 
found for the VPM in each of these four cases, we can combine the results into a single Jones 
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matrix. 


J(d) 


J(d) 


cos (j) — siiup \ / 1 0 \ ( d) cos(f> r Eh (0,d) sirup \ 
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This Jones matrix maps the incident electric field ^ E™ E te j to the reflected field 

^ E™' E TEl 1 . The Stokes’ transfer function of the system can then be determined as a 
transformation of the density matrix describing the transverse correlations in the system, 

15' = J ] 15J. ( 9 ) 


4. Circuit Parameter Values 

The circuit topology described above is physically motivated; however, circuit parameters 
for the grid are required to analyze the response of the VPM. Many treatments found in 
the literature for wire grid polarizers consider a simple analytically-treatable limiting case 
[8,12,27,28]. These models assume no diffraction, A < 2 g, little or no azimuthal dependence 
of current density on the wires, a <C A, and that the grid filling fraction is small, 2a/ g < 
1 / 27 T. However, in the millimeter and submillimeter, tension exists between realizable grid 
structures and those that strictly meet these conditions. A survey of practical grid polarizers 
indicates that filing fractions in the range 0.3 < 2a/ g < 0.8 are available and achieve 
high polarization isolation. Such grids are outside of the range of applicability of the models 
referred to above. Green’s function [28] and lattice-sum [29,30] methods enable more accurate 
treatments without a need for these approximations. Here we employ a hybrid numerical 
approach that allows us to analyze the behavior of wire grid polarizer used for this work. In 
this section, we first revisit the classic analytical treatment of polarizing grids and motivate 
the need for a more rigorous treatment. 


4- A. Sparse Grid Approximation 

As an example of commonly-employed analytical methods, we explicitly give the expressions 
for the circuit elements from Marcuvitz [12]; however, this is representative of the approxima- 
tions used in quasioptical treatments of grids. In the following, all reactances are normalized 
with respect to the impedance of free space. Given a ware radius a and a separation or 
“pitch”, g, the circuit elements for the polarization parallel to the wares (the “inductive” 
mode; Marcuvitz §5.21) are given by 
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2 


( 10 ) 



The series for the inductive reactance in Equation 11 is rapidly convergent. In the limit of 
applicability, this series can be approximated by ~ 0.6 • (r;/ A) 2 and is a small perturbation. 
As a result, the factor in brackets, which we identify with the inductance, is effectively 
frequency-independent. 

For the polarization having the electric field perpendicular to the wires (the “capacitive” 
mode; Marcuvitz §5.13), the reactive circuit elements are given for a symmetrically-centered 
capacitive post [12]: 


X c = 

Xcc = 



1 


4 



where A\ and A 2 for normal incidence are defined as 


( 12 ) 

(13) 



The circuit parameters given above are plotted as functions of the geometric filling fraction, 
2a,/ g, in Figure 4. Leading wavelength-dependent terms have been factored out in order to 
isolate circuit elements as a function of grid geometry alone. This separation is anticipated 
in the quasistatic limit. The triangles at the bottom of the figure indicate typical values 
of grids employed in the submillimeter and millimeter, and the vertical gray line marks the 
position at which 2a/ g = 1/n. We note that the expression for the shunt inductance becomes 
negative at this point. This is a manifestly unphvsical result, as the wires do not actually 
become capacitive for the parallel polarization. Thus, a different approach is required in this 
regime, as the limit of applicability indicated by Marcuvitz, 2a/ g < 1/2tt, is violated. 
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Fig. 4. The Marcuvitz circuit parameters are plotted as functions of the ge- 
ometric filling fraction, 2a/ g. The vertical gray line is shown at 2njg — l/n. 
The shaded region indicates the region in which the Marcuvitz approxima- 
tion ceases to hold [12]. Triangles indicate typical values of 2 a/g for grids 
constructed for submillimeter and millimeter use. 
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4.B. Numerical Simulations of Grid Performance 

In order to understand the inductance and loss of wire grids, we perform numerical simula- 
tions as a function of filling fraction over the range 0.02 < 2a/ g < 1. The numerical grid sim- 
ulations were carried out with CST (Computer Simulation Technology) MICROSTRIPES™ 
Time-Domain Transmission-Line Matrix (TLM) solver and Ansoft HESS (High Frequency 
Structural Simulator). See Figure 5 for the model geometry for each polarization and the 
boundary conditions employed to define the structure. We consider loss from the structure 
in the limit where the field penetration depth is small compared to the wire diameter and 
the incident wavelength is greater than two times the wire pitch. In this quasi-static limit, 
the higher order Floquet harmonics or grating lobe responses are absent, and the field con- 
figuration is a perturbation on the lossless case. 

Wheeler’s incremental inductance rule [22] can be applied to compute the resistive losses 
from the wires in the grid. This approach implicitly assumes that the current distribution on 
the conductors does not vary appreciably over distances comparable to the thickness of the 
wire and is a convenient parameterization for data reduction. This condition is satisfied for 
a wire having a well-developed skin effect, S <C a, in the single mode limit. The anticipated 
form can be expressed as a product of two separable functions: the first dependent on the grid 
filling factor, 2a/ g, and the second a function of skin depth over the wavelength, 5/ A. Guided 
by this physical insight, we proceed to model the grid in the single mode approximation to 
further explore the scaling properties of the grid inductance and loss. 





Fig. 5. The unit cells used for the simulation in both the inductive (left) and 
capacitive (right) modes are shown. An infinite structure is represented by the 
symmetry of the boundary conditions. 

The power reflection coefficient, R , from the impedance presented by the shunt reactance 
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of the grid, Zl = Rl + iX L , is: 
R 


u 1 2 _ 

*->11 — 
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(l + 2R L y + (2X L f 


where all impedances are normalized to Z 0 . The resulting reflection phase is 
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(16) 


(17) 


For a near optimal grid geometry, the impedance contrast with respect to free space, 2 Z L <C 
1, and the resultant reflection phase is approximately 7r. Similarly, the power transmission 
can be expressed as 


2Z l 2 (2R L f + {2X L f 
1 + 2 Z L ~ (1 + 2R l ) 2 + ( 2 x L y ' 


(18) 


and the power absorption, A = 1 — R — T = 1 — |5n| 2 — l^il 2 , is computed from conservation 
of energy. Here, Sij are the simulated complex scattering parameters for the structure. A 
corresponding set of expressions can be written for the capacitive case. We note the ratio of 
the absorptanee over the reflectance can be expressed as 
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4 R s 
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(19) 


The bulk surface resistivity of the metal is Rs — 1 /aS, r] is the grid loss efficiency, and A is 
the observation wavelength. The field penetration depth, 5 = yj 2/ /icru) , is computed from 
the metal’s permeability, //,, electrical conductivity, a, and the observation frequency, u>. ft 
is informative to consider the dimensionless ratio that occurs in the circuit model above in 
terms of this parameterization: 


4 Rs 4\]ujfi/2a 4 ttS 



( 20 ) 


The expression on the right is related to the Hagen- Ruben emissivity formula [23], and can 
be derived from Fresnel coefficients and identified as the emissivity from a bulk metal at 
normal incidence. This approximation is valid in the limit of low emissivity. Equation 20 
links the circuit theory to the interaction between the materials and the electromagnetic 
fields. We see that the loss efficiency, rj — (A/R)(4 tt6/X)~ 1 , is the loss of the grid relative to 
a flat sheet made of the same material. We compute the scattering parameters as a function 
of grid filling factor and frequency and derive equivalent circuit parameters. In Figure 6 we 
show the grid loss factor, the grid shunt reactance, and the grid reflection phase appropriately 
normalized to remove the wavelength dependence as a function of grid filling factor. This 
scaling is anticipated in the quasi-static limit for the circuit elements. 
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Fig. 6. (TOP) The grid loss factor is plotted as a function of geometric 
filling factor for a wire grid. Both HFSS and Microstripes simulations are 
shown for the cases of the electric field perpendicular to and parallel to 
the grid wires. The filling factors for typical polarizing grids are plotted for 
comparison. (MIDDLE) The reactance for the inductive mode is shown as 
a function of the geometric filling factor. A useful interpolation function is 
X n = 0.51 — 1.19\n(2na/g) + 0.53(ln(27ra/f/)) 2 + 0.11(ln(27ra/</)) 3 . For com- 
parison, the Marcuvitz model is shown as well as the measured inductance for 
the Hertz grid. (BOTTOM) The normalized reflection phase of the inductive 
mode is plotted as a function of geometric filling factor for the HFSS simula- 
tion. In this case, we find (p n = 1 — \/3.545 {2a/ g). In each of the panels, the 
condition 2 n/g — l/V, is denoted by a dashed vertical line. Triangles indicate 
typically manufactured grid geometries. 
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For the perpendicular illumination case, we find a grid loss factor, g± = 2 a/g. For the 
parallel illumination case, some care is required- the conductive wires effectively define a 
finite thickness aperture. For this configuration, the influence of the waveguide cutoff has a 
non-negligible influence on the grid’s absorption properties, and we find its influence results 
in a grid loss factor of = 7/tanh(77r • 2 a/g) for this polarization. In referring to Figure 6 we 

note that this differs from previous grid loss factors reported in the literature. For example, 
for infinitely thin strips, Ulrich [24] finds g\\ = g/2a. For small filling factors, this expression 
and the one derived above have the same slope; however, in the limit 2 a/g -7 1, the influence 
of the exclusion of the fields by the grid wires has a pronounced influence on the response. 
This behavior is also observed in our numerical simulations. We note that within the context 
of the planar waveguide model [22], by symmetry, the field configuration for this polarization 
is physically equivalent to that of a microstrip via [25] . 

The numerical simulations enable us to calculate the inductance in the range of interest 
to compute the response of the grid. The inductance is plotted as a function of 2 a/g in 
the middle panel of Figure 6 and compared to the Marcuvitz model. Our numerical study 
predicts the shunt inductance asymptotically approaches zero for large values of 2 a/g and 
does not become negative as observed in the Marcuvitz approximation outside the stated 
range of applicability [12], 

Numerous authors have indicated that the “optimal” grid spacing satisfies g = 2?ra. A 
degree of caution needs to be exercised regarding this point. Namely, the derived inductance 
expressions are of logarithmic accuracy, and higher order terms are required in this regime, 
e.g., see Landau §33 [31], and [32]. This is demonstrated in the middle panel of Figure 6 
where we have simulated the structure and extracted the reactance. We caution in pursuing 
an expansion in 2 a/g as a means of analysis, many terms are required before S/X and «/ A 
are relevant. Our numerical approach avoids this issue. 

In realizing an optimal wire grid polarizer, one must consider the effect on both the 
parallel and perpendicular polarizations. For the parallel polarization, larger 2 a/g leads to 
high reflectance since the shunt inductance tends to zero in this limit. For the perpendicular 
polarization, larger 2a / g leads to a decrease in the magnitude of the capacitance between 
wires and undesirably increases the reflection of this mode. Grids having g ~ 2na perform 
well because the wire spacing is small enough that the inductive mode is highly reflective, 
but large enough such that the capacitive mode is highly transmissive. 

5. 350 /im VPM Tests 

5. A. Experimental Setup 

Preliminary characterization of the Hertz VPM instrument has been reported [2], Our setup 
is similar to that used previously and is schematically depicted in Figure 7. A chopped black- 
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body is polarized by a wire grid and focused onto the Hertz polarimeter that simultaneously 
measures horizontal and vertical polarizations. The VPM under test is placed near a pupil 
in a collimated beam within the fore optics that couple the source to the Hertz cryostat. For 
a given grid-mirror separation, the polarization is modulated by rotating a cold half-wave 
plate (HWP) in discrete steps. The aperture of the blackbody is sized to roughly fill a single 
Hertz detector. Inside of the Hertz cryostat, the polarization is separated into two linear 
orthogonal components by a wire grid analyzer. Both the transmitted (T) and reflected (R) 
components are detected by bolometers. 



Hertz Cryostat ^ 

"T" Bolometer 
"R" Bolometer 


-0 


"tea 


* 

Analyzer I 


Polarizer Rotating 
□ Chopper 



Folding 

Flat 


Modulator Throw 


VPM 

(wires parallel to page) 


Fig. 7. The main elements of the experimental setup for the 350 //m test 
are shown [2], Radiation from a blackbody source is polarized by a wire grid 
polarizer having wires oriented at an angle of 45° with respect to the plane 
of the page. A chopper modulates the intensity of the signal. The radiation 
is collimated prior to being processed by the VPM. Upon exiting the VPM, 
the radiation is relayed to the Hertz cryostat. Inside, the radiation passes 
through the half-wave plate and the bandpass filter before being diplexed into 
two orthogonal linear polarizations. A bolometer detects the signal in each 
polarization. 


The polarization transfer function of the VPM is measured as follows. At various grid- 
mirror separations, the output linear polarization is measured using the Hertz polarimeter. 
The normalized Stokes parameters are then reported as a function of the grid-mirror sepa- 
ration. 



5.B. Data Analysis 


At each grid-mirror separation, the basic data analysis pipeline is similar to that described 
in [16]. However, we describe the process in some detail since there are some deviations 
that are associated with the difference between laboratory measurements and those done 
astronomically. At a single HWP position, 16 “frames” are recorded. Each frame consists of 
two co-added "chops” or on-off differences. In our laboratory setup, the chopper phase was 
observed to drift relative to the commanded TTL signal from the data system. To correct 
for this, the in-phase data and the quadrature data were both used and the signals for a 
single frame in the R and T arrays are 


riioej) = 

U(aj) = 


in phase 
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(21) 

(22) 


Here, oq is the value of the HWP angle at the jth HWP position. This technique mitigates 
the problems associated with phase drift in the chopper. 

The frames recorded for each HWP are averaged separately for each of the R and T signals, 


%)4E r ‘W < 23 ) 

V i=l 
1 JV 
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i=i 

Here, N = 16 is the number of data frames taken at each HWP position. This process is 
repeated for each of six HWP positions (nominally, 0°, 30°, 60°, 90°, 120°, and 150°). The 
relative gain between the R and T bolometers is calculated as the ratio of the sum of all the 
data in the R bolometer and that of the T bolometer: 


E, «(*i) 

E jTiajY 


(25) 


The polarization signal is then calculated using both the R and T bolometers. This combi- 
nation aids in the removal of common-mode noise in the system, 


s ( n j 


jjM - /Hw) 
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(26) 


For each grid-mirror separation, the polarization signal in Equation 26 is fit to the function 


S(a) = Qcos(4o) + U sin(4a) + V\ cos(2a) + l 2 sin( 2 a). (27) 
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The parameters Q and U are the linear Stokes parameters. We use the following convention 
for normalized Stokes parameters: q = Qj I , a = U / 1 , and v = V/ 1 . The V Stokes parameter 
is related to V\ and V 2 via |V| = \J (Vj 2 + V 2 2 ). The system was found to be dominated by 
drifts on the time scale associated with the HWP cycle. To account for this, error estimates 
are determined by goodness of fit of the data to Equation 27. 

The inclusion of the circular polarization in the half-wave-plate merit function was found to 
be required; however, this is not unexpected. Half-wave plates are perfect half-wave devices at 
a single frequency that is dependent on the birefringence of the crystal (An) and its thickness. 
Deviation from perfect matching between the HWP center frequency and that of the band- 
defining filter as well as the details of the input spectrum can both cause some leakage 
between circular and linear polarizations. In this case, we have found that the mixing occurs 
at the level of ~ 20% (i.e. for an input circular polarization of unity, one would measure 
v = 0.20). At this level, the reduction of efficiency of the HWP as a linear polarization 
modulator is small- the residual linear polarization efficiency is only reduced by ~ 2%. This 
is consistent with previous characterization of the Hertz polarimeter [17,18]. Because of the 
modulation topology of the VPM. the HWP was exposed to much larger amounts of v than 
in an astronomical instrumentation setting, thus emphasizing this otherwise small effect. 

Initially we observed that the q Stokes parameter exhibited a modulation as a function 
of grid-mirror separation that had a similar modulation pattern to that of Stokes v. This 
effect was observed at the 15% level. By rotating the Hertz cryostat by 90° and remeasuring, 
this effect was found to switch signs, indicating some birefringence in either the lenses or 
the pressure window in the cryostat. Since the axis of such a birefringent element could a 
priori be oriented in any orientation, a similar effect may be contaminating the u Stokes 
parameter signal. To mitigate this possibility, the reported values for the VPM transfer 
function were evaluated by averaging each data point with an equivalent data point taken 
in a configuration where the cryostat was rotated by 90°. The resulting values for q and u 
are shown as a function of grid-mirror separation in Figure 8A. 

We have applied the circuit model to data taken using the Hertz VPM instrument with a 
single VPM having wires nearly parallel to the plane of incidence. In this case, it is necessary 
to average over all relevant frequencies: 

/ ‘OO 

D(u) j\v)ij) n {v)dv , (28) 

-OO 

where i> n (v) is the normalized instrument bandpass (approximately uq = 353 gm, Xv/v = 
0 . 10 ). 

The q and u data are fit to the VPM model using the downhill simplex method [19]. Ten 
fit parameters were varied in this minimization: the angle of incidence, an offset in the grid- 
mirror separation, a global polarization efficiency, a HWP phase offset, Xlc, Xi, Rl, Xc, 
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Xcc i and the bandwidth stretch (k). The bandwidth stretch takes into account deviation 
from parallelism between the grid and mirror in the VPM. Such an effect would cause the 
bandpass to appear slightly wider than that defined by the bandpass filter. 

The best fit parameters for our circuit, model arc given in Table 1. Where applicable, 
these values are compared with the values estimated from the system setup, the Marcuvitz 
circuit values, and/or the circuit values derived from our numerical simulations for the grids 
(a — 12.5 /tm and g = 67.5 /mi.). The resulting models for q and u are plotted along with 
the data in Figure 6A. The reduced \ 2 value for the joint (q,u) fit is 2.1. Upon removal of 
five of the (q, u) pairs where either q or u has a significant deviation from the model, the ,\' 2 
is reduced to 1.2. From this we conclude that the circuit model captures the basic features 
of the polarization transfer function, however, the statistical errors are underestimated in a 
subset of the measurements. Such an effect can occur as a result of time- variable instabilities 
in the system that are insufficiently mitigated by the techniques described earlier in this 
section. Based upon the observed performance of the system a phase drift in the chopper 
wheel for a subset of the data frames is the most probable cause. 

The resistance in the inductive mode, Hi, is higher than expected from the formulation 
given above. The expected value is based on the conductivity of gold, the coatings for the 
wires. This is an indication that the gold coatings are compromised; the fitted resistance is 
consistent with the base metal (tungsten, p = 1.89 x 10 7 fl • m). The fitted value of Xl is 
within 30% of that derived from the numerical simulation described above. We have plotted 
this value shown in the middle panel of Figure 6. The model presented here depends on the 
ability to identify a unit cell for the grid; however, variations in wire spacing violate the 
symmetry used to calculate the properties of the grid using this technique. Examination of 
the grid using a microscope confirms this concern. 

In Figure 8B, the features in Stokes q at grid-mirror separations of approximately 200 
and 400 pm are due to trapped resonances for the inductive mode. We note in passing that 
these features can be significantly reduced by modifying the VPM configuration. In place 
of the mirror, one could incorporate a second grid with wires oriented perpendicular to the 
first followed by beam dump. Ideally, the brightness temperature of the beam dump should 
match that of the source over the beam on the sky. This architectural change would reduce 
common-mode systematic effects between orthogonally-polarized detectors. 

5 .( 7 . VPM Modulation Efficiency 

The modulation efficiency of the Hertz polarimeter was experimentally observed to be ~ 95%. 
This measurement was made with the source grid in place but in the absence of the VPM 
polarizer. This value is consistent with previous measurements of the efficiency of the Hertz 
polarimeter [2,18]. For the data reported in Table 1 and Figure 8, the reported efficiency, 
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Grid-Mirror Separation (am) 

Fig, 8. (A) The q and u values were measured as a function of grid-mirror 
separation at 350 fin 1 . (B) Measurements of the Hertz VPM at 3 mm are also 
shown. The transmission line model is in close agreement with the infinite 
wavelength approximation at this frequency; however, the difference between 
the two is measurable. 
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Parameter 

Estimate 

Fit 

Units 

Angle of Incidence 

20 

22.83 

[degrees] 

Grid- mirror Offset 

0 

-2.01 

[/an] 

Polarization Efficiency 

100 

98.4 

[%] 

II WP Phase Offset, 

0 

0.0125 

[radians] 

K 

1 

1.08 

[-] 

Hl 

9.9 x 10~ 4t 

3.4 x nr 3 

H 

A 'l 

0.031 t 

0.040 

[-] 

X lc 

-0.26* 

-0.27 

[-] 

X c 

- 

-0.43 

[-] 

Xcc 

- 

-265 

M 


Table 1. Fit results for the 350 fi in measurements. Fro in §4, *From Marcuvitz) 


98.4%, is that of the VPM alone. The lower modulation efficiency previously reported, 85% 
[2], was found to be reproducible and induced by beam misalignments at the location of 
the VPM. The estimate reported here for the modulation efficiency of the VPM is a lower 
limit and could arise from one of several effects. First, the circuit model assumes uniform 
wire spacing and size. As such, the model represents the average circuit parameters over the 
illuminated part of the grid and does not account for variations in parasitic coupling between 
the wires that are important at high frequencies. Treatment of this effect would require 
modifications to the simple circuit model considered here. Alternatively, and more likely, the 
non-unity efficiency factor could be an artifact of the system such as a residual misalignment 
of the beam on the VPM or residual uncertainty in the Hertz efficiency calibration. 

6. 3 mm VPM Tests 

We have also measured the VPM at a wavelength of 3 mm that corresponds to gj A = 0.022. 
We measured the response from 75-110 GHz; however, we used the 91-100 GHz portion of 
the measured spectrum to synthesize a 10% bandwidth to match the 350 /mi measurement. 
The setup for this test is shown in Figure 9. An HP 8510 vector network analyzer was used 
to measure the polarization state of the radiation as it is processed by the VPM. Vertically- 
polarized radiation is emitted from port 1 and launched from waveguide to free space via 
a feed horn. It is then reflected by a polarizing grid having wires oriented in the vertical 
direction. This grid helped to further define the polarization state of the radiation and limits 
the formation of trapped reflection resonances in the system. The unused polarization is 
terminated by a free-space beam dump. An ellipsoidal mirror is used to form a beam waist 
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on the VPM. This ensures that the phase front of the radiation is approximately flat at the 
polarization modulation point. .Just as for the 350 /mi setup, the VPM wires were oriented 
such that their projected angle is 45° with respect to the incoming polarization state. A 
second, identical ellipsoidal mirror refocuses the beam on a second feedhorn that is identical 
to the first. A folding flat is used in this final free-space path to symmetrize the optics. The 
orthomode transducers (OMTs) [20,21] ports corresponding to the unused polarizations are 
terminated by precision waveguide loads. This mitigates the influence of polarization mode 
conversion that can result in trapped resonances in the beam waveguide structure. 

As the grid-mirror separation is varied, the VPM modulates the incoming (vertical) signal 
that we define as Stokes U in order to maintain consistency with the 350 /mi data. The 
modulated power is detected by port 2 which is sensitive to vertically-polarized radiation. 
The 800-point spectrum of 5 12 over the entirety of W-band (75 to 115 GHz) is recorded at 
every grid-mirror separation. Once this is finished, a 90° twist is added to the waveguide in 
port 2 to switch the port’s sensitivity from vertical to horizontal polarization. The grid-mirror 
separation of the VPM is then adjusted in an identical way as before. 

6. A. Data Analaysis 

For comparison with the 350 /mi data, we wish to determine the response of the VPM over 
a 10% bandwidth fractional bandwidth. To do this, we calculate the integrated signals H(z ) 
and V(z): 


H(d) = J2 |Sf 2 (d ,%)| 2 


(29) 

(30) 


Here, S^(d, v) and S\ 2 {d, v) are the relevant scattering parameters for the case of horizontal 
and vertical sensitivity for port 2, respectively. 

The next step is to calculate the relative gain of the two channels. This can be done by 
noting that the single- frequency response for the two quantities are, 


H ( rf ) = \[?h + Uh cos 4>(d) + V H sin 0(d)] (31) 

V(d) = - I v — Uy cos (j){d) — Vy sin 0(d)]. (32) 


Here, I h/v ■ Uh/v , and Vh/v are Stokes parameters. 

For each of H (d) and V(d). we perform a linear least-squares fit. Since the single-frequency 
will not fit our broadband response precisely, we do not use this to directly calculate Stokes 
parameters. Rather, we calculate a relative gain based on the unmodulated component 
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Beam 

Dump 


T.T ~~ t Modulator 

T Throw 

VPM 

(wires at 45 deg.) 


Fig. 9. The setup for the 3 mm VPM transmission test is shown. Radiation is 
emitted from port 1 of the HP 8510 with a vertical polarization. The radiation 
is then reflected off of the grid and collimated by an ellipsoidal mirror. At this 
point, the radiation is reflected off of the VPM. An identical ellipsoidal mirror 
follows, and a folding flat directs the beam into a second feed that is attached 
to port 2 of the HP 8510. A 90° twist is added to the waveguide to change 
the sensitivity from the vertically-polarized state to the horizontally-polarized 
state. 
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(Stokes I) from each fit, 


(33) 


From this, Stokes u can be calculated: 


/ 


hi 

V 


= m-jm 

[) n(d) + fv(dy 


(34) 


The results of this analysis are shown in Figure 8B. In addition, we have plotted a VPM 
model using reactances appropriately scaled from those found in the 350 gm fit, using the 
leading order wavelength dependence expected from the Marcuvitz circuit elements. We 
found that to fit the data required an incident angle of 19° and an offset in the grid-mirror 
separation of -30 /mi. No other parameters in the model have been adjusted. For comparison, 
we have plotted the normalized Stokes parameter, u , in the infinite-wavelength limit. In this 
limit, the impedance contrast is high for the inductive mode and the dependence oir the 
circuit parameters is minimal. 


7. Discussion 

Figure 10 summarizes the VPM transfer function for a single frequency for each of the 
models in this work as well as the geometric limit. To emphasize the differences in shape 
of the curves, for each case we have shifted the curve to align the first minimum to lie at 
zero delay. For A m> a. the grid circuit is firmly in the single-mode limit. In this case, for 
plane wave illumination, Equation 2 is a good approximation for the VPM phase delay. As 
the wire diameter becomes a finite fraction of a wavelength, the impedance contrast of the 
inductive circuit is reduced. The polarization response remains a sinusoidal function of the 
phase delay; however, the VPM reflection phase is dependent upon the details of the grid 
geometry and impedance contrast. 


8. Summary 

We have applied a transmission line model to the VPM used for the Hertz/ VPM submil- 
limeter polarimeter and have obtained estimates for key circuit parameters using numerical 
simulations. The simulations represent an improvement over analytical approximations in 
estimating the inductance of wire grid polarizers having 2wa ~ g, such as that used in this 
work. 

At 350 //m the VPM grids have wires that are ~0.07A in diameter, corresponding to a 
free-space phase shift across the width of a wire that exceeds 10°. We have shown that in 
this regime, it is not sufficient to identify the VPM phase delay with that derived from 
the path difference induced by the grid-mirror separation. The VPM transfer function is 
modified from the infinite wavelength case, both in shape and in position. This is due to a 
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combination of the reactances imposed by the thickness and spacing of the wires and the 
interaction of these reactance with the variable position of the mirror. At 3 mm, where the 
impedance contrast for the inductive mode is higher, our measurements demonstrate much 
closer agreement to Equation 2. 




— 



g/X 

0.19 

0.022 

0 

a/X 

0.036 

0.0042 

0 

X (gm) 

350 

3,000 

OO 


Fig. 10. The modeled VPM reflection phase delay is shown for single- 
frequencies at 350 and 3000 gm for the grid models obtained above. In the 
limit g / A <C 1, a sinusoidal form for Stokes U is observed. As this condition is 
relaxed, the VPM reflection phase delay differs from the free-space grid-mirror 
delay. 
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